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Abstract 



We show that when Bjorken xb is large, the leading order twist-four contri- 
butions to the structure functions can be expressed in terms of the derivatives 
of the normal twist-two parton distributions. Our analytical results are not 
only consistent with the parameterized power corrections extracted from the 
data, but also predict the flavor dependence of the power corrections. We also 
discuss the impact of our analytical results on the extraction and the flavor 
separation of the parton distributions at large x. 
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I. INTRODUCTION 



Perturbative Quantum Chromodynamics (QCD) has been very successful in interpreting 
and predicting scattering cross sections at large momentum transfer. For hadronic collisions, 
a number of sets of reliable parton distributions are available for perturbative calculations 
of potential experimental signals |l|]. However, recent data on inclusive high Pt jets from 
Fermilab @] and high Q 2 events at HERA || remind us that there are still large uncertainties 
in the parton distributions at large x. The precise knowledge of parton distributions at large 
x is very important in search for signals of new physics at very high energies, which might 
be reached only by the partons with large momentum fraction x at the colliders. 

According to the QCD perturbation theory, parton distributions at large x and high Q 2 
can be derived from the parton distributions at large x and lower Q 2 by solving the DGLAP 
evolution equations. However, data from most physical observables are only sensitive to 
the parton distributions in the intermediate or small x region, because of the phase space 
integration and the fact that the parton distributions decrease very fast as x increases. It 
is believed that structure functions measured in deeply inelastic scattering (DIS) at large 
Bjorken xb and low Q 2 might be the best source to get the direct information on the 
parton distributions at large x. But, at large xb and low Q 2 , the structure functions receive 
significant power corrections [||, and these power corrections confuse the determination of 
the parton distributions. 

In general, there are two types of power corrections, known as the kinematic and the 
dynamical power corrections. The kinematic power corrections is due to the fact that the 
nucleon has a finite mass, rriN- Such power corrections have been systematically discussed 
in Ref. ||. On the other hand, the dynamical power corrections are the results of the 
interactions between quarks and gluons The dynamical power corrections are in 

general proportional to the multiparton correlation functions or the matrix elements of high 
twist operators, which are different from the normal parton distributions. Because of the 
large number of unknown multiparton correlation functions, the actual size of the dynamical 
power corrections to the structure functions have not been well understood theoretically. In 
this paper, we derive analytically the leading dynamical power corrections to the structure 
functions, and provide the numerical estimates of the power corrections. Since we will not 
discuss the kinematic power corrections in this paper, without any confusion, the words 
"power corrections" in the rest of this paper indicate the "dynamical power corrections" . 

In a spin-averaged lepton-hadron deeply inelastic scattering, l{k) + h{p) — * l{k') + X, the 
hadronic tensor W^^XByQ 2 ) can be decomposed as || 



where q = k — k' is the four-momentum of the virtual photon, q 2 = —Q 2 , and Fi(xs, Q 2 ) with 
i — 1,2 are the structure functions. At the large enough Q 2 , the structure functions can be 
factorized into a convolution of the short- distance coefficient functions and nonperturbative 
matrix elements of multiparton fields f|, 




(1) 



2 



F t (x B ,Q 2 ) = Y,C^ m \x B ,x 1 ,x 2 ,..,Qy^)®T^ m \x 1 ,x 2 ,.., f i 2 ) , (2) 

where Xi with % = 1, ... are parton momentum fractions, fj 2 is the factorization scale and is 
often chosen to be Q 2 in DIS, C^' m ' are the perturbatively calculable coefficient functions, 
and TV' m > are the corresponding nonperturbative matrix elements. In Eq. J2 m runs 
from to infinity, and J2j sums over all possible types of nonperturbative matrix elements 
at a given power of m. The term with m = is often known as the twist-2 or leading 
power contributions, and the terms with m ^ are the high twist contributions or power 
corrections. The factorization presented in Eq. (^) is useful only when the terms of high 
power corrections can be neglected. When the physical scale, Q 2 , is large enough, the leading 
power contributions are often sufficient, and the structure functions can be expressed in terms 
of a few nonperturbative matrix elements, T^'°\x, Q 2 ). These matrix elements T^'°\x, Q 2 ) 
with j — q,q and g are known as the normal quark, antiquark, and gluon distributions, and 
have been well-determined for a wide range of x and Q 2 HJ. 

Working in the moment space of the structure functions, -Fj(n, Q 2 ) = f dxx n Fi(x, Q 2 ), 
De Riijula, et al. || demonstrated that when the moment n is large enough, 



F 2 (n,Q 2 )~£ 



Q 2 



T m (n,Q 2 ), (3) 



with /ig a non-perturbative mass scale, and concluded that there exist a region of n ~ Q 2 / fi 2 , 
where the power corrections are as important as the leading contributions. This conclusion 
is equivalent to say that in the x^-space, when x B —>■ 1, high order terms in the expansion in 
Eq. (H) become important even if Q 2 is large. Physically, when x B — ► 1, the invariant mass of 
the photon-hadron system in DIS approaches to the resonance region, and the perturbative 
expansion in Eq. (0) is then invalid. 

In order to avoid the resonances, we have to require the invariant mass of the photon- 



hadron system ||,[T0| 



W 2 = (p + q) 2 = m 2 N + Q 2 > 2 GeV 2 . (4) 

x B 

On the other hand, in order to extract the information on parton distributions at large x B 
and low Q 2 , we have to work in a region where W 2 is small and not too far away from the 
resonances, and will have to deal with the power corrections. The purpose of this paper is 
to get a good estimate of the size of the leading power corrections and to discuss its impact 
on the extraction of the reliable twist-2 parton distributions at large x. 

By fitting DIS data in large x B region, it has been found flJTT| that the 1/Q 2 dynamical 
power corrections to the structure functions can be parameterized in terms of the normal 
parton distributions, 

F 2 (x B , Q 2 ) « (l + y&^J F 2 (x B , Q 2 ) LT , (5) 

where F 2 (x B , Q 2 )lt includes all leading twist (or leading power) contributions as well as the 
target mass corrections to the full structure function F 2 (x B ,Q 2 ). The h(x) in Eq. (|^) is a 
phenomenological fitting function, and is parameterized as [111]] 
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Ii(xb) = a 



x 



B 



1 - X B 



— c 



(6) 



where a, b and c are constant fitting parameters. Because of the 1/(1 — Xb) dependence in 
Eq. (||), it is clear that the dynamical power corrections in Eq. (|5|) become more important 
when Xb increases. It has been argued that the simple parameterization in Eq. (|6]) can 
provide a very good fit to the existing DIS data on proton as well as deuteron targets 
(with slightly different fitting parameters), and demonstrated that the power corrections are 
important for extracting the correct parton distributions, such as d/u ratio, in large x region 

IB 

On the other hand, a complete expression of the 1/Q 2 power corrections to the structure 
functions at the leading order of a s were derived more than fifteen years ago [ 12 , 13 1 . As 



expected from Eq. (|j), the 1/Q 2 corrections are proportional to the matrix elements of four- 
parton operators, which are very different from the operators defining the normal parton 
distributions. The success of the phenomenological parameterization in Eq. (||) raises a 
question: under what conditions, one can derive the approximate expression in Eq. (||), or 
a similar expression, from the complete 1/Q 2 power corrections at the leading order of a s , 
which were derived in Ref. []12| , ^3| . 

In this paper, we show that when xb is large, the 1/Q 2 power corrections to the structure 
functions at the leading order of a s can be approximated as 



F 2 (xb,Q 2 ) ~ e q X B 

q,q 



1 + 



Q 2 



(D 2 T ) 



d 



-x B 



dx 



B . 



+ 4 (m 2 T ) 



q{x B ,Q 2 



(7) 



where m 2 N = p 2 is the nucleon mass, (D 2 -) and (m\) are the averaged values of D\ and m^, 
which will be defined in the next section. Numerically, we demonstrate that our analytical 
expressions for the leading power corrections in Eq. (^) have the same x^-dependence as 
the phenomenological parameterization given in Eqs. (||) and @, and therefore, we expect 
Eq. (0) to fit the data as well. Instead of three unknown parameters, a, b and c in Eq. (|), our 
analytical result in Eq. (|7|) have only two parameters, and both of them have clear physical 
interpretations. Furthermore, because of the scaling violation of the quark distributions, 
the derivative in Eq. (0) gives the natural explanation for the flavor dependence (or target 
dependence) of the fitting parameters in Eq. (||) [|TTf| . In addition, the derivative provides 
some extra logarithmic Q 2 -dependence for the power corrections. 

Our result in Eq. (|7|) demonstrates that although the power corrections to the structure 
functions in DIS are in general very complicated, the leading contributions in a s at large x b 
are approximately given by the derivatives of the normal quark distributions. Such a simple 
analytical expression will be very useful for the QCD global analysis of extracting reliable 
leading twist parton distributions. 

The rest of our paper is organized as follows. For the completeness, we briefly derive the 
leading twist-4 contributions to the structure functions in Sec. f|. In Sec. [TTI|, we identify 
the necessary approximations, and derive the main result of this paper, Eq. (]?]), from the 
complete leading power corrections obtained in Sec. [TT]. Finally, in Sec. [TV] , we numerically 
compare our analytical expression in Eq. (^) with the phenomenological parameterization in 
Eqs. (|^) and (|B|), and estimate the values of the parameters (D 2 -,) and {m 2 -) in Eq. (|7p. We 
also discuss the possible impact of our analytical results of the leading power corrections to 
the extraction and the flavor separation of parton distributions at large x. 
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II. POWER CORRECTIONS IN DIS AT LEADING ORDER IN a s 



The complete 1/Q 2 power corrections to the structure functions at the leading order of 
a s can be separated into two groups: one is proportional to the two-quark-two-gluon matrix 
elements (or correlation functions), and the other is proportional to the four-quark matrix 



elements [12,13]. Because of Fermi statistics and the fact that sea quark distributions are 



smaller than the gluon distributions, contributions from the four-quark subprocesses in large 
xb region should be smaller than that from the quark-gluon subprocesses. In this section, 
we neglect the four-quark contributions, and drive the leading order 1/Q 2 contributions to 
the DIS structure functions from the quark-gluon subprocess, from which the main result of 
this paper, Eq. (|7|), will be deduced in the next section. 

For the defmiteness of our derivations, we choose the following photon-hadron frame 

Q 2 

= —xbP ■ nn^ H , 

XbP ■ n 

pt 1 = p ■ nff , (8) 

where we have neglected the hadron mass for the dynamical power corrections. In Eq. 
the four-vector, n M and n M are defined as 



nl 1 



IV 



(n+,n",n T ) = (1,0, T ) 



n" = (0,l,0r) , (9) 
n 2 = 0, n 2 = 0, and n • n — 1 . 

In terms of the n^ 1 and n M , we can reexpress the hadronic tensor in Eq. ([[]) as 

W^(x B ,Q 2 ) = ef F L (x B ,Q 2 ) + eg F 1 (x B ,Q 2 ) , (10) 

where the longitudinal and transverse tensors are defined as 



e L U = 77^? [x B p-nn^ + — — nA [x B p-nn v + — — 

2Q Z \ zxsp-n J \ ixsp-n 

eg = n^n" + n v - g" v = dT , (11) 

with ex eg = 1/4, ex v eg = 2, and ereg = 0. The longitudinal structure function, 
F L (x B ,Q 2 ) in Eq. flT^), is defined as usual, F L = F 2 /x B - 2F 1 . From Eq. ([10]), the DIS 
structure functions can be extracted from the hadronic tensor as 

F 1 (x B ,Q 2 ) = ^egW llu (x B ,Q 2 ) 

F L {x B ,Q 2 ) =4ef W^xb.Q 2 ) . (12) 

As derived in Ref. [|13[], the 1/Q 2 power corrections to the structure functions from the 
quark-gluon subprocesses can be factorized, as shown in Fig. [T], 



W^(x B ,Q 2 ) 



l/Q 2 Q 



^- J dxdx 1 dx 2 Tr Cg^(x B ,x,x 1 ,x 2 )T a(3 (x,x 1 ,x 2 ,Q 2 ) , (13) 
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where "Tr" is the spinor trace for fermions. The partonic part, C^Jxb, x, x±, x 2 ) in above 
equation, is given by the top parts of the Feynman diagrams in Fig. |. The blobs 111 
Fig. |T]a are given by the diagrams in Fig. 0a, and the blobs in Fig. |I|b and Fig. [TJc are given 
by the diagrams in Fig. |||b. In Fig. |2|, the fermion lines with a short bar indicates that 
the propagator includes only the short- distance contact term, as defined in Ref. The 
quark-gluon correlation functions, T a/3 (x, x\, X2, Q 2 ) in Eq. (0), is expressed in terms of the 
four-parton matrix element, 



T^(x,x 1 ,x 2 ,Q 2 ) = J 



p + dy p + dy l p + dy 2 



i{x-x 1 )p+y 1 -i(x-x 2 )p+y 2 



27T 2tt 2tt 



(14) 



where D a (y) = d a + igT A B (y) is the covariant derivative, which differs by a factor of 



% from the definition used in Refs. [ 12"lJi~3" |. In Eq. (0), Dtp 
line-integrals between the fields are omitted ]13| . 



(9 



a/3 



n 



a n p )D Pl and the 



To complete the factorization between the partonic parts and the twist-4 matrix elements, 
we need to separate the spinor trace, as well as the Lorentz indices linking between the C^i 
and T af3 in Eq. ([13]). Such factorization can be achieved by decomposing the T a/3 as 



T afi = \l-pd^T+ l - 1 -p lb (ie^) f + 



(15) 



where T and T are scalar functions of x, Xi,x 2 and Q 2 , and ". . . " represents terms that 
are further power suppressed. In Eq. fll5]), the symmetric tensor d a/3 is defined in Eq. ( O) 
and the antisymmetric tensor e 



a/3 



n,,n 



„. From Eqs. dl2|) and fll3|) , we can derive the 
1/Q 2 contributions to the structure functions at the leading order in a s by calculating the 
Feynman diagrams in Fig. [|. Substituting Eq. (|i"5|) into Eq. (|I3~D , and contracting e T u with 
the hadronic tensor W^ v (x B , Q 2 )\i/q 2 , we obtain in the light-cone gauge (n ■ A = 0), 



e^W^{x B) Q 2 
et?W$(x B ,Q 2 



1/Q 2 






2 2x B 


1/Q 2 


6q Q 2 




2 2x B 


1/Q 2 


6q Q 2 



dx 1 dx 2 6{Xl - XB) T q (x 2 ,x 1 ,Q 2 ) 



dxi dx 2 



x 2 - 
5(x 2 



x B 
- x B ) 



X\ — X B 



T g (x 2 ,xi,Q 2 ) 



(16a) 
(16b) 
(16c) 



where e q is the quark(or antiquark)'s fractional charge of flavor q, and the superscripts (a), (b) 
and (c) correspond to the contributions from diagrams Figs. 0a, 0b and 0c, respectively. The 
quark-gluon correlation function T q (x 2 , x\, Q 2 ) is defined as 



T q (x 2 ,x 1 ,Q 2 



dy P + dy 1 



e iX!p+y e i(x 2 -x 1 )p+y 1 



2tt 2tt 

x (p\4>Ml + (d af sD a (y^D^y^) ^ q {y~)\p) . 



(17) 



Substituting Eq. (|16D into Eq. (| 

2 



F!(X B ,Q' 



1/Q 2 



we obtain 
<5(x 2 



dxi dx 2 



x B ) -8(x 1 - x B ) 



x 2 - Xi 



T q {x 2 ,x u Q 2 ) . (18) 
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Similarly, by contracting to the hadronic tensor W^(x B ,Q 2 ), we obtain 

ef W${x B ,Q 2 ) =e\ [ dx 8{x — x B ) T q (x, Q 2 ) , 
^ i /o 2 CJ J 

el u WV(x B ,Q 2 ) 



e?W$(x B ,Q 2 ) 





= e 2 — 


i/Q 2 


q Q 2 




= , 


i/Q 2 






= , 


i/Q 2 





where the quark-gluon correlation function T q (x,Q 2 ) is defined as 



T q (x,Q 2 



1 r dy 



4 J 2tt 

Combining Eqs. (|i~2"D and (|l~9|) , we have 



e ixp+y ~ (p|^(0) 7a ^(0) 7 + 7/3^(l/-)^(l/-)|p) • 



Fl(x b , Q 2 ) 



VQ 2 g 



9 Q 2 



T q (x B ,Q 2 ) . 



(19a) 
(19b) 
(19c) 



(20) 



(21) 



From Eqs. (|i"8|) and (|2l|) , we obtain the leading order 1/Q 2 contributions from quark- 
gluon subprocess to the structure function F 2 (x B ,Q 2 ) as 



F 2 {x B ,Q 2 ) 



= Yl e l x b\ 2x b / dx x dx 2 
i/Q 2 Q T 



<5(x 2 - x B ) - 5(x 1 - x B ) 



-8 Jdx5(x-x B )T q (x,Q 



x 2 - x 1 

2\ \ 



T^x^x^Q 2 ) 
(22) 



Our results given in Eqs. (|T8"D, (pl|) and (0) are consistent with those in Refs. [0,0, after 
taking into account the difference in the definition of the covariant derivative (a factor of 
V) and an extra "-1" in the definition of T q (x, Q 2 ) in Eq. (|0|). 

III. LEADING POWER CORRECTIONS AT LARGE X B 

In this section, we identify the approximations that enable us to derive Eq. (|7]) from 
Eq. (|22"D, which represents the complete 1/Q 2 contributions to F 2 (x B ,Q 2 ) from the leading 
order quark-gluon subprocesses. 

From the phenomenological fitting of the dynamical power corrections |H|] , 1/Q 2 con- 
tributions to the structure functions are most important in large x B region, because of the 
1/(1 — x B ) factor in Eq. (^|). From the complete leading order 1/Q 2 contributions in Eq. (^2|), 
it is nontrivial to conclude that the 1/(1 — x B ) behavior exists. However, under some simple 
approximations, we demonstrate below that Eq. (^) has indeed the 1/(1 — x B ) behavior 
when x B is large. 

The second term in Eq. (EE3) is simple and directly proportional to the quark-gluon 



correlation function T q (x B ,Q 1 ) in Eq. 



To see if the T q (x B ,Q ) has the 1/(1 — x B ) 



behavior, we simplify the operator expression of the T q (x B ,Q 1 ) by using the equation of 
motion, , y a D a (y)ip q (y) = 0, and we have 
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T q (x B ,Q 2 



1 

~4 
(x B ) 



dy~ 



2n 



(p|^(0)7- J D + (0)7 + T - J D + ( 2 /-)V g (2/-)b) 



p + dy 
2n 



ix B p^y 



<pfe(o) 



p+ 7 - 



)Ip) 



(23) 



where the light-cone gauge n • A = and D% = (g a/3 — n a n^)Dp were used. On the other 
hand, the normal twist-2 quark distribution can be expressed in the light-cone gauge as 



q(x B ,Q 2 ) 



dy- 



Jx B p+y 



27T 



m 



N 



p + dy 
2vr 



tx B p^ y 



V 7 



i> q {y )\p) 



(24) 



where m 2 N = 2p + p~ is the nucleon mass square. In order to relate the correlation function 
T q (x B , Q 2 ) to the quark distribution q(x B , Q 2 ), we estimate the difference between the matrix 
elements (jp\^ q {0)\p + 'y~]if>q{y~)\p) in Eq. (|3[) and (j>\ip q {ti)[p~l + }4> q (y~)\p) in Eq. as 
follows. For free field operator ijjqS on free light quark states, the difference between these 
two matrix elements is the difference between [p + k~] and \p~~ k + ] with k the light quark 
momentum. In our frame, we can choose 



p = m 2 N /{2p + ) and k = m 2 ^/{2k + ) 



and k~ 

,2 , 7,2 



x B p 



(25) 



with the transverse mass of the light quark to be m\ = m q + k^. If we assume the free field 
relation can be generalized to the full fields, we obtain 



p + dy 
2n ~ 

m 2 N (x B f 



jxbv y 



(p|^(o) 



PI 



)Ip) 



p + dy 
2tt 



(p|^(0) 



P 7 



)Ip) • 



(26) 



Substituting Eq. ( p6| ) into Eqs. (|23| ) and (p4[), we derive 



Q 2 ) « - {m 2 T ) q(x B , Q 2 ) 



(27) 



We emphasize here that the importance of Eq. ( f2~7j ) is not how accurate the approximation is, 
rather it indicates that the second term in Eq. (|22|) does not have the 1 / (1— x B ) enhancement. 

The first term in Eq. (f22j) depends on the correlation function T q (x2,x\,Q 2 ), which can 
be sketched in Fig. [| and be reexpressed as 



T q (x2,xi,Q 2 



dy_ 

2tt 

p + dyl 



2tt 



e ^v + y~ (p|^(0) 7 + 



i> q {y )Ip) 



(28) 



Because of the factor [5{x2 — x B ) — 5{x\ — x B )]/(x2 — x\) in Eq. fl22|) , x\ ~ X2 ~ x B in 
T q (x2,xi,Q 2 ) when is large, and the effective "gluons" in Fig. [5] are soft. We then 
assume that the effective "gluon" operator d a pD a (yi)D f3 (yi) in Eq. (^8|) is not sensitive to 
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the position variable y 1 , and further approximate this operator by an averaged expectation 
value as 



d aP D a (y^(y^ « (d a pD a (0)D^(0)) = (D 2 T ) 



(29) 



More discussion on this assumption will be given in the next section. With the approximation 
in Eq. (E9), we can rewrite the correlation function T q (x 2 , Xi, Q 2 ) as 



T q (x 2 ,xi,Q 2 



dy 
2tt 



~(D 2 T )8{x 2 -x x ) 



i> q {y )\v) 



(D 2 T )5(x 2 -x 1 ) 



(30) 



where q(xi,Q 2 ) is the twist-2 quark distribution defined in Eq. (PH) . Using Eq. (|30|) , we 
derive 



J dxi dx 2 



5(x 2 - x B ) - 5(xi - x B ) 
x 2 - Xi 



T q (x 2 ,xi,Q 2 



dxs ) 

Substituting Eqs. (g7|) and (gl]) into Eq. 
the structure function F 2 (xb,Q 2 ) as 



If 1 

Z JO 

q(x B ,Q 2 ) . 



d 



(31) 



F 2 (x B ,Q 2 



i/Q 2 



1 



I>9 X ^ 



, we obtain the leading power corrections to 
d 



dx 



+ 4 (m 2 T ) 



B . 



(32) 



Combining this leading power corrections and the leading order twist-2 contributions to 
F 2 (xb,Q 2 ), we obtain our main result, Eq. (|7|). 

Although Eq. ([32] ) looks very different from the phenomenological parameterization in 
Eq. (H), we show below that Eqs. (|32|) and (||) have the same analytical behavior as xb — » 1- 
When xb is large, the quark distribution have the following behavior 

q(x B ,Q 2 )^{l-x B y^ asz B -l, (33) 

with a positive and real parameter n q . In general, n q (Q 2 ) has the logarithmic Q 2 -dependence 
due to the scaling violation of the parton distributions. For the valence quark distributions, 
n q ~ 3. Consequently, we have 



d 



- x B - 



-q{XBj 



X R I 

n q(Q 2 ) ^ ~ q( x BiQ 2 ) + terms without 



' dx B iy yv ~ v/ i_ i B m -v / 1 — Xs 

Substituting Eq. (|34"D into Eq. (0), naturally, we derive the leading 1/(1 — x B ) behavior, 

1 



(34) 



F 2 (x B} Q 2 



i/Q 2 



x B q{x B ,Q 2 



(Di) ,^ 2 . x B 

' - n q (Q ) h terms without 



Q 1 



X B 



(35) 

which is consistent with the phenomenological parameterization in Eq. (§) that was intro- 
duced in Refs. after fitting the existing DIS data. 
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IV. NUMERICAL RESULTS AND CONCLUSIONS 



In this section, we demonstrate numerically that as a function of xb, our analytical 
expression in Eq. ([?]) and the phenomenological parameterization defined in Eqs. (|5|) and 
(0) have almost identical behavior, even though they look very different. 



In the following numerical estimates, we use CTEQ3L parton distributions of Ref. [Tj] for 
the quark distributions in the proton. In addition, in order to avoid the resonance region and 
to test the perturbative power expansion of the structure functions, we keep W 2 > 2 GeV 2 
in our calculations. 

In Fig. [|, we plot the structure function F 2 (xb, Q 2 ) as a function of ig for Q 2 = 4 and 
9 GeV 2 . The solid line is the lowest order contribution, F 2 °(xb,Q 2 ) = J2q£qXBq(xB,Q 2 )- 
The dotted and dashed lines are given by the expressions in Eqs. (|5|) and (|7|), respectively, 
with F 2 °(xb, Q 2 ) for the F 2 (xb, Q 2 )lt in Eq. (|5|). For the fitting parameters, a, b and c in 
the phenomenological parameterization h(x), we use the following values a = 0.50, b = 3.2 
and c = 0.11, which are given in Ref. ||11|| . For our analytical expression in Eq. (|7|), there 



are two unknown parameters (-D 2 -) and (m 2 ^). We choose (-D 2 -) = 0.1 GeV 2 to be consistent 
with the overall size of the power corrections. The transverse mass square (m 2 ^) for the light 
quarks should be of the order of the intrinsic transverse momentum square, and therefore, 
(my) > Aq CD . The actual size of (m\) should depend on the collision energy W 2 due to the 



soft gluon shower from the incoming quark line ||15|| . Within our approximation, different 



choices for the value of (m 2 ^) represents a constant shift of the size of the leading power 
corrections in (F 2 — F 2 °) / F 2 ° . In the rest of our discussions, we choose (m^) to be zero for 
a close exam of the derivative term. In this paper, we are more interested in understanding 
the functional dependence of the power corrections than trying to extract the best values 
of the (D 2 -) and (m^) by fitting any data. We defer the detailed global QCD analysis for 
extracting the parton distributions from the data to another publication. 

Because of the steep falling parton distributions when xb increases, we plot the relative 
power corrections to the structure function, (F 2 — F 2 °) / F 2 ° , as a function of xb in Fig. [5|. 
As shown in Fig. |5], the leading order power corrections are very significant and can be 
as large as 30% in large i s region. In addition, with only one parameter (D 2 -,) (because 
(m 2 ^) is set to be zero here), our analytical result in Eq. (|7|) is remarkably consistent with 
the parameterization extracted from the data. The difference between the solid and the 
dashed curves in Fig. |5| cannot be separated by the existing data. Such consistency on 
x^-dependence is already a nontrivial test for the QCD perturbation theory. 

As shown in last section, the leading power corrections at large xb are given by the 
derivatives of the normal quark distributions. It is the derivative that makes the power 
corrections different for quark distributions of different flavors. We introduce 



,2, 1 



Aq{x B ,Q 2 ) = — 



(D 2 T ) (- XB J-)+A(m 2 T ) 
\ dx B J 



q(x B ,Q 2 ), (36) 



for quark flavor q. In Figs. |6| and [7], we plot the relative power corrections Aq/q as a function 
of xb for the "up" and "down" quark, respectively. Because the "down" quark distribution 
falls off faster than the "up" quark distribution in the CTEQ3L parton distributions, the 
"down" quark obtains a larger power correction than what the "up" quark gets. Such 
difference is clearly evident in Figs. ^| and [7|, and is very significant when xb is large. The 
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power corrections derived here also receive some logarithmic Q 2 -dependence through the 
derivative of the parton distributions, due to the scaling violation of the twist-2 parton 
distributions. 

As demonstrated above, the power corrections to the structure functions are not only 
significant at large xb, but also sensitive to the quark flavors. It is such flavor difference 
of the power corrections that makes the extraction and the flavor separation of the parton 
distributions nontrivial, in particular, in large xb region. Recently, a lot of attention have 
been devoted to the ratio of d/u ||. We defer the detailed global QCD analysis for the 
extraction and the flavor separation of the parton distributions to another publication. 

In summary, we derive the analytical expressions for the leading power corrections to 
the DIS structure function, given in Eq. ([32]). We made two main assumptions in our 
derivation, which are introduced in Eqs. fl26|) and (^), respectively. Since the 1/(1 — xb) 
enhancement of the power corrections is from the term proportional to the quark-gluon 
correlation T q (x2, X\, Q 2 ), the assumption in Eq. ( P9"| ) is more crucial than the other. In 
general, it may be possible that the (Dj) is not a constant independent of the position (or 
equivalently the value of Xb)- However, as long as any part of the (D 2 -) is not proportional 
to (1 — xb), we will have the 1/(1 — xb) enhancement. We believe that at large xb, at 
least the "soft" gluon part of the averaged value of the covariant derivative square should 
not be sensitive the location inside the nucleon, even though the averaged value of the 
transverse momentum square may be proportional to (1 — Xb)- Therefore, we conclude that 
the leading order power corrections to the structure function F2(xb,Q 2 ) should have the 
1/(1 — xb) enhancement, and our analytical result, Eq. (]?]) with two fitting parameters 
(Dj) and (rrij>), should predicts the correct leading x^-dependence of the power corrections 
at large xb- Because it has only two unknown parameters, our result has a good predicting 
power and can be tested at CEBAF as well as future fixed target programs at Fermilab. 
It is also useful for the QCD global analysis for extracting parton distributions in large xb 
region. In addition, our approach can be applied to the Drell-Yan process to study the 
leading power corrections there fll6| . 
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(a) 




(b) 



Fig. 2 



FIG. 2. Feynman diagrams for the blobs in Fig. |T[ The wave lines and curly lines represent 
the virtual photons and gluons, respectively. The solid lines are for the quarks, and the line with 



a short bar represents the contact term of the propagator [13] 




Fig. 3 



FIG. 3. Diagram represents the quark-gluon correlation function T q (x2,x\,Q 2 
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Fig. 6 

FIG. 6. Relative power corrections to the "up" quark distribution at Q 2 = 4 GeV 2 and 
Q 2 = 9 GeV 2 , respectively. 
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Fig. 7 

FIG. 7. Relative power corrections to the "down" quark distribution at Q 2 = 4 GeV 2 and 
Q 2 = 9 GeV 2 , respectively. 
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